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Abstract. We provide a new approach to studying the Dirichlet-Neumann map for 
Laplace's equation on a convex polygon using Fokas' unified method for boundary value 
problems. By exploiting the complex analytic structure inherent in the unified method, 
we provide new proofs of classical results using mainly complex analytic techniques. 
The analysis takes place in a Banach space of complex valued, analytic functions and 
the methodology is based on classical results from complex analysis. Our approach 
gives way to new numerical treatments of the underlying boundary value problem and 
the associated Dirichlet-Neumann map. Using these new results we provide a family 
of well-posed weak problems associated with the Dirichlet-Neumann map, and prove 
relevant coercivity estimates so that standard techniques can be applied. 
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1. Introduction 

In the last twelve years there has been a rapid development of the so-called Fokas 
method for boundary value problems [13J. This method was initially developed for the 
study of boundary value problems associated with integrable nonlinear PDEs [9l [HJ [12] . 
Remarkably, the method has proved an extremely powerful tool in the study of linear 
boundary value problems [H [10] and has offered new avenues of pursuit in the numerical 
study of such problems [HI [151 ESI EE] . 

The Fokas method can be informally considered as the Fourier analogue of the 
classical boundary integral methods [171 [19]. In the latter case, the analysis is done in 
the "physical" space - i.e. the space associated with the domain on which the boundary 
value problem is posed. In the Fokas approach the analysis is done in "spectral" space 
and the classical boundary integral equations are replaced with the global relation. Like 
the usual boundary integral equations from potential theory, the global relation gives a 
relationship between the known boundary data and the unknown boundary values for 
a given PDE. However, the form of the integral equation is very different: the global 
relation has meromorphic dependence on a spectral parameter which plays the analogue 
of the wave number in Fourier analysis. The Dirichlet and Neumann terms that arise 
in boundary integral equations have spectral analogues that arise in the global relation. 
The fact that the spectral boundary data has meromorphic dependence on the spectral 
parameter allows for use of powerful tools from complex analysis and this salient feature 
is one of the core reasons behind the success of the Fokas method. 

Another major ingredient in the Fokas approach is the use of novel integral 
representations. By formulating and solving Riemann-Hilbert problems related to the 
underlying boundary value problem, Fokas has developed a means of representing the 
solution to many important problems. These integral representations are intimately 
related to the classical fundamental principle of Ehrenpreis and Palamodov [HI 120] . 
A direct consequence of this highly abstract result is that any solution to a constant 
coefficient PDE on a convex domain can be written as a superposition of exponential 
solutions. More concretely, if P = P(—id/dxi, . . . , —id/dx^) is a constant coefficient 
differential operator and Pu = in a convex domain in H N , then the fundamental 
principle states that u can be represented in the form 

u{x) = [ c{x, X)e iX - x dfi(X) 
J z P 

where Z P = {A G C N : -P(A) = 0} and c(x, ■) is a polynomial in x such that c(x, X)e lX ' x 
for A G Zp solves Pu = 0. This theorem is highly abstract and non-constructive. 
The result states that there exists a measure d/i supported on Zp and a collection 
of exponential solutions c(A, x)e lX ' x such that the representation holds. Fokas' novel 
integral representations are realisations of this abstract result - providing an explicit 
expression that is precisely a superposition of exponential solutions. These integral 
representations must, of course, contain terms relating to the known boundary data. 
The novelty lies in the fact that these terms arise in the same form as they do in the 
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global relation. This means there should be no need to go back and fourth between 
physical space and spectral space - it is sufficient to understand the global relation in 
spectral space and give the solution in terms of integrals of spectral functions by means 
of the novel integral representation. 

The implementation of the Fokas method has been largely formal in nature. One 
usually works on the assumption that a solution to the underlying problem exists and 
aims to construct it by analysing the global relation. This assumption can then be 
checked a posteriori. Recently more rigorous results have been obtained [21 El H]. In 
this paper we continue to address rigorous aspects of the Fokas approach. 

We concern ourselves with the particular case of Laplace's equation in a convex 
polygon. For Dirichlet data with square integrable tangential derivatives along the 
edges of Q, we prove the following results (more precise versions of which can be found 
in Theorems El [3] and their corollaries) 

(I) The global relation defines a continuous linear map from the spectral Dirichlet 

Data to the spectral Neumann data. 
(II) The resulting physical Neumann data is square integrable along each edge. 

(Ill) The global relation gives rise to an infinite family of well-posed weak problems 
that are easily approximated using standard Galerkin techniques. 

As a corollary to statements (I) — (III) we get a new proof of existence for the classical 
Dirichlet problem in a convex polygon for boundary data in H 1 (dQ) [27]. The statement 
in (II) is in accordance with the well known regularity of the Steklov-Poincare operator 
on Lipschitz domains [15] . In proving (III) we provide a basis from which the previous 
numerical studies [HI [T5| [251 126] can be made rigorous. The methods presented here 
can easily be adapted to deal with the Neumann boundary value problem. In this case 
similar statements (I) — (III) hold, but with the Dirichlet data being determined modulo 
constants. 

The approach we use is far removed from the classical methods boundary integral 
and more modern pseudodifferential methods. The analysis takes place Banach space of 
complex analytic functions and the main ingredients in our proofs are complex analytic 
in nature, utilising the classical theorems of Liouville, Montel, Phragmen and Lindelof. 

2. The Global Relation for Convex Polygons 

We work on a polygon Q C R 2 ~ C with vertices {£j}™ =1 and sides r« = (zi, Zi+i) with 
z n+ \ = z\. We write on = arg(^ + i — Zi) for the angle the side r, makes with the positive 
real axis and Ay = oij — aj. Also set 2a = for the length of the side IY 

We are given real valued boundary data fi G if 1 (Pj) for i = l,...,n, meaning 
that fi and its first tangential derivative are square integrable along the edge Tj. We are 
interested in the Dirichlet-Neumann map associated with the classical Dirichlet problem 



Aq = in fi, 
q = fi on T { for i = 1, . . . 



n. 



(la) 
(lb) 
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That is to say, we want to reconstruct the unknown Neumann boundary values from 
the known Dirichlet boundary data. 

In practical applications the more physically relevant quantity in (1) is the gradient 
field, Vq = (dq/dx, dq/dy), or equivalently the complex derivative dq/dz with z = x+iy. 
In [8] it was shown that any solution to (II ah has an integral representation 



dz 2tt J u 



where the spectral functions {pi(A)}™ =1 are defined by 



A(A) 



r, 



dz 1 



dz' 



(2) 



(3) 



and the {^j^Li are rays in the complex plane orientated out towards infinity with 
arg(A|f i ) = —on- The spectral functions satisfy the global relation 



!>(*) = o. 



(4) 



8=1 



Note that the spectral functions ([3]) contain information about the known boundary 
data and the unknown boundary values since dq/dz' involves derivatives in the 

tangential and normal directions along IV We can interpret (j2J) as a formal solution 
to (1) if we assume the spectral functions satisfy the global relation Our aim is 
to solve Pj) for the unknown parts of the spectral functions which contain information 
about the unknown normal derivatives. 

It was shown in [15] (c.f. [3j Hj) that on the assumption that the global relation 
(jl]) is satisfied, the integral representation (j2]) provides a solution to the boundary value 
problem (1). This is important from both the theoretical and practical point of view - it 
means that solving the global relation for the unknown parts of the spectral functions is 
equivalent to solving the boundary value problem (1). By describing the global relation 
(TjJ as a map between function spaces for the spectral functions, we are able to provide 
existence, uniqueness and stability results for the solution to the global relation. Perhaps 
more importantly, this gives a means for the practical numerical solution to the global 
relation for the unknown parts of the spectral functions. 

It will be convenient to have a local description of the edges IY Let us introduce 
the local parametrisations ipi : [— <Ti,(Ti] — > IV with 



Mr) 



1 

2ct~ 



(cr, + r)z i+1 + (a; - r)zi 



= rrii + re ' 



where m ; = \{z,i + is the mid-point of the side IV For a function / : Tj — > C 
we write its pullback by ipi by ip*(f){r) = f{jp{j)). Using this notation the spectral 
functions are written 

dq 



A(A) 



-iXrrii 



We note that 



Oq 



dq . dq 



dz 



re 



dr. 
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where d/dt and d/dn denote the tangential and outward normal derivatives along Tj. 
Setting (p\ = ip*(dq/dt) and (pf = i/j*(dq/dn) we can write the spectral functions as 

„— iXmi 



^(e^A) + i^ n (e^A) 



where we have defined the Fourier transform 



J= : ^ i-)> tpi(X) 



-i\r 



<Pi(r) dr. 



The global relation then takes the form 



8=1 



-iArrii 



^ n (e iai A)-i^(e iai A) 







(5) 



which holds for all A G C. It will be convenient to have a more symmetric form of the 
global relation. To this end, fix some i G {1, . . . ,n}. Multiply (jSJ) by e lXmi and replace 
A with \e~ loli . We find 

#(A) - i#(A)] + ^ e ie " iaiA ^-^) [^(e-^A) - ^(e-^A)] = 0, 

for each 1 < % < n. Set $ n = ■ ■ ■ ,<PnY> ^ = {0u ■ ■ ■ > ^Vf an< i introduce the 
operator T = I + K, where I is the identity and K is the linear operator defined by 

$i(A) H- (K$)i(A) = ^e ie ~ iai(mi -^ )A $ J (e- iA «A) ) 1 < i < n. (6) 

Then the global relation can be written succinctly as 

T($ n - i$*) = 0, A € C. (7) 

Each of these n equations are equivalent to the original global relation (jSJ). The vectors 
$*(A) and $ n (A) contain the spectral boundary data, which in this case is just the 
Fourier transform of the original functions. In what follows we characterise spectral 
Dirichlet-Neumann map $* h- > $ n defined by (J7J). 



3. The Real and Complex Paley- Wiener Spaces 

Here we discuss the relevant function spaces that will be used in the sequel and cement 
some of our notation. 

The global relation for Laplace's equation has been given in (J7J). The components 
{</3*}™ =1 of the known vector $* are are related to the derivatives of the Dirichlet data 
fi G H 1 (r i ), and we have <p\ G L 2 [— <7j, crj for i — 1, . . . , n. The global relation contains 
the Fourier transform of this data. It is natural then to work with the classical Paley- 
Wiener spaces 

which contain the Fourier transforms of square integrable functions whose support is 
contained in the interval [— <7i,<7i]. The classical Paley- Wiener theorem states 

PW ai = {/: C C entire, / \f(x)\ 2 dx < oo, |/(A)| < e e am+e) Ve > o), 
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i.e. the space PW ai consists of entire functions of exponential type <Tj whose restrictions 
to the real axis are square integrable. Paley- Wiener functions satisfy the important 
pointwise inequality 

l/(*b)l £ ll/llae" 1 * 1 , fePW*. 

We will refer to this as the Paley- Wiener inequality. For a standard treatment of the 
Paley- Wiener spaces we refer the reader to j5j HHJ ED EI] . 

It will be convenient to work with slightly modified versions of the Paley- Wiener 
space. In the general setting of the Paley- Wiener theorem the space L 2 [— <jj, <Tj] refers 
to complex valued, square integrable functions. However, our data will be manifestly 
real. Hence forth we shall use 0i, <xj to denote the space of real valued, square 
integrable functions. We will then work on 

where the subscript "sym" refers to symmetric. The reason for this is highlighted in the 
following simple lemma. 

Lemma 1. The space PW° y l m is a closed subspace of PW a ' whose members obey the 
symmetry condition 



f(X) = /*(A) = /(-A) 

for all A G C. 

Proof. That PW^ m is closed in PW Ul is clear. In addition, if / G FL\[— <7j, then 
the symmetry condition is satisfied. Conversely, if / G PW ai and obeys the symmetry 
condition then 



= /(A) - /(-A) = ^ e-^{g{r) - ^(r))dr 



for some g G L 2 [— <jj\. But then the Fourier inversion theorem implies that 
9 e L%[-o-i, (J;], so / G TX^J-o-j, <Ti]. □ 

Remark 1. It is clear that the classical Paley- Wiener space can be decomposed as 
PW° l = PW a% ffi PW a% 

sym ^ ' ' asym 

where the latter space consists of Fourier transforms of imaginary valued, square 
integrable functions. The characterisation of this space is given by an anti-symmetry 
condition, where an extra minus sign appears. 

It is well-known [24, Ch. 6] that PW Ci is a closed subspace of L 2 (R). It follows 
that PW ai is a Hilbert space when equipped with inner product 



(/i,/ 2 )= / fx(x)f 2 (x)dx. 



The same is true of PW^ m with this inner-product. This means that one can treat 
PW^ m as a closed subspace of £ 2 (R) whose elements have an analytic extension to the 
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entire complex plane which is of exponential type <7j and obeys the necessary symmetry 
condition. 

We introduce the function space X = PW ai x • • • x PW Un with norm 

V2 / „ \ 1/2 



( " /-CO \ / 



where here and throughout || ■ H2 denotes the usual L 2 norm on the real line. Using the 
decomposition PW ai = PW^ m © PW^ ym we may write X as 

X -A S y m © -A aS y m . 

It will be convenient to regard X sym as the "real part" of X, while X asym is the "imaginary 
part". Both are Banach spaces in their own right when equipped with the norm || • 
Also set Y = L 2 (Rr) xn with norm 

V2 / n \ 1/2 



|2 
12,- 



where here and throughout || • H2,- denotes the usual L 2 norm on the negative real axis. 
It is straightforward to show that || • \\x and || ■ ||y are equivalent norms on Agy m and 
Aasym owning to the symmetry and anti-symmetry properties of the elements of the 
respective spaces. We note the isomorphism X ~ L 2 (dQ). Each of X sym , X asym and Y 
are Hilbert spaces when equipped with the appropriate inner product, but we shall only 
need their Banach space structure. 

We will often refer to the Fourier transform of an £ 2 (R) function, and this is to 
be understood in the limit-in-the-mean sense [22]. If T : U — > V is a continuous linear 
map between normed spaces we write T e C{U,V). A norm-bounded subset S of a 
normed space U is one in which there is some constant C such that \\u\\ < C for each 

u e S c u. 



4. Some Functional- Analytic Results 

Here we prove some functional-analytic results for the Paley- Wiener spaces which will 
prove useful for the purposes of studying the spectral Dirichlet-Neumann map. 

Throughout this section X will denote an arbitrary measure space with positive 
measure /i. We use L P (X, (i) with p e [1, 00) to denote the Banach space of (equivalence 
classes of) complex valued measurable functions on X with norm 

f ^ (X' /|Pd/i 

The following theorem will be of particular importance when studying the continuity of 
the Dirichlet-Neumann map, but seems of interest in its own right. 

Theorem 1. Let T : PW a — > L p (X,fi) be a continuous linear operator. Suppose also 
that T is also continuous with respect to the topology of point-wise convergence, i.e. if 
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{fn}n>i is a sequence in PW a and f n — )■ / pointwise, then Tf n — > Tf pointwise in 
L p (X,fi). Then T has closed range. 

The proof of requires the following "pseudo-compactness" lemma which which will 
be of use throughout the paper. 

Lemma 2. Any norm-bounded sequence in PW U contains a subsequence that converges 
pointwise and locally uniformly to an element of PW U that obeys the same norm bound. 

Proof. Let {f n }n>i be norm bounded in PW° . Then using the inequality Paley- Wiener 
inequality |/0o)l < e CT ' 2 ° ||/|| 2 we deduce that the sequence {/„} n >i is locally uniformly 
bounded, i.e. for each compact K C C we have sup^ |/ n (;z)| < Ck for some constant 
Ck- By Montel's theorem, we can extract a convergent subsequence that converges 
pointwise, locally uniformly to some analytic function. Let {f nk }k>i be this subsequence 
so that for each compact K 

lim sup|/ njk (z) - f(z) \ = 

k— >oo k 

for some analytic function /. We claim that / e X with the same norm bound. First, 
note that / is certainly of exponential type a. Indeed, if we fix R > and choose k 
sufficiently large so that \f(z) — f nk (z)\ < 1/Rior \z\ < R we have 

SUp |/(»| < SUp |/0) - fn k (z)\ + SUp \fn k (z)\ 
\z\<R \z\<R \z\<R 

<^ + \\fn k he^ 
< e alzl 



where the final constant is independent of R. The norm bound is an immediate 
consequence of Fatou's lemma. □ 

Remark 2. One might examine the statement of this lemma and speciously reason that 
the Paley-Wiener spaces are Montel spaces, i.e a uniformly bounded subset of PW a 
contains a convergent subsequence. This would be false because PW a is a Banach 
space, and since the unit ball is not compact in an infinite dimensional Banach space 
it cannot possibly be a Montel space. What we have shown is the following: given 
a bounded sequence {/ m } m >i in PW a , one can extract a subsequence {/ mfc }fc>i that 
converges pointwise, locally uniformly to some / e PW a . However, we have not shown 
that f mk — > f in PW a , i.e. ||/ — /mjh — > 0, and in general this will not be the case. 
Indeed, if we consider the standard basis functions for PW a 

/o r sin(a"A — nm) _ /o r (— l) m sin(crA) 
\n (aX — nm) \n (aX — irm) 

then it is clear that ||/ m ||2 = 1 but f m — > pointwise and locally uniformly. 

Proof of Theorem U\ It is enough to prove that T maps norm-bounded, closed sets in 
PW a to closed sets in L p (X,[i) [TJ p. 79]. Let Xg be a closed, norm-bounded subset 
of PW a with || /|| 2 < 5 for each / E Xg. Set g n = Tf n for some sequence {/ n } n >i 
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in X$ and suppose g n — > g in L p (X,n). Then it is well-known that one can extract a 
subsequence {g nk }k>i that converges pointwise almost /i-everywhere to g [231 Th. 3.12]. 
Hence 

g(x) = lim (Tf nk ) (x) almost /i-everywhere in X. 

The sequence {/n fc }fc>i is & norm-bounded sequence with each ||/n fc ||2 — 3- By Lemma 1, 
we can extract a subsequence {f nkl }i>i that converges pointwise (and locally uniformly) 
to some / e Xg. Using the pointwise continuity of T we deduce 



g(x) = lim [Tf n 



T ( lim f n ) (x) = (Tf)(x) almost /i-everywhere in X. 



So g = Tf in L p (X,[i) for some / G X$. Hence the image each norm bounded closed 
set in PW a is closed in L p (X,[i), and we deduce that T must have closed range. □ 

5. The Spectral Dirichlet-Neumann Map 

Recall the global relation ([7]) is 

T($ n - 1$*) = 0, A e C 

where T = I + K, with K defined in ([6]). We will need some properties of this operator. 

Lemma 3. We have T e C(X, Y) and T is also continuous with respect to the topology 
of pointwise convergence. 

Proof. That T respects pointwise convergence is obvious. To prove the relevant estimate 
for the first claim we use J 7-1 ^ = I on £ 2 (R) to write 

2tt ./ ,, 

Using the convexity of the domain Q one sees that for j ^ i, i ± 1 
e < arg ( e" 1Ql (re 1Qj + m j - m*) ) < n - e 



for some e > 0. After an application of Cauchy-Schwarz and Parseval's theorem we get 
estimates of the form 



g ie ,a *(mi-mj)A<£^ e -iAy^ 

So for A < we have 



< e Asine ||$ 



J II 2 - 



|(T$)j(A)| < |$i(A)| + 



+ e~ |A|sin l$|| x . 



The || ■ || 2,- norm of the first and last terms are clearly dominated by ||$||x; so we need 
only look at the remaining terms. Using the representation OH]), these terms are, after an 
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appropriate change of variables, the Laplace transform (along a ray) of an L 2 function 
supported on [0, 2a j] for j — i ± 1. They have the generic form 

e^ A / ' e-W^hjWdr, (A < 0) 
Jo 

for appropriate hj and Wj with Rewj > 0. It is well-known [16] that the Laplace 
transform defines a bounded linear map from L 2 (R + ) to L 2 (R + ). We deduce that the 
|| • || 2,- norm of these terms are bounded by constant multiples of ||$j||2 for j — i ± 1. 
Using Minkowski's inequality we find ||(T$)j|| 2i - II^IU- By applying these estimates 
to each of the n components of T$, we deduce that T G £(X, Y). □ 

Remark 3. Obviously T G £(X sym , Y) and £(X asym , Y) also. 

Remark 4. In the context of the previous lemma, we should interpret the global relation 
(J7J) as describing a map from X to Y . In this case we should only really consider A G R . 
However, all the terms appearing in (JTj) are entire functions so we can make a unique 
analytic extension to the entire complex plane. We play fast and loose in this regard, 
making no distinction between T$ defined on R~ and its analytic extension, say (T$) ext , 
defined on the the entire complex plane with (T$) ext = T$ on R~. 

In this current setting it is clear that the spectral Dirichlet-Neumann map $* i— > $ n 
is determined by the null space of the operator T G C(X, Y). We write N(T) for the 
null space. We arrive at the following problem: 

Given $* G X sym find $ G N(J) with Re$ = $*. 

The Neumann data is then $ n = lm$. We note the analogy between the operators 
T -<->■ 8, the latter being the <9-derivative which annihilates complex analytic functions. 
The functions $* and $ n playing the roles of the real and imaginary parts of the analytic 
function. For each $ G X sym we set 

DN($) = {$' G X sym : $ + i$' G N(J)}. 

This set will prove useful. 

Lemma 4. For each $ G X sym the set DN($) is a singleton. 

Proof. First we prove that DN($) contains no more than one element, and this is 
equivalent to showing that T is injective on X sym . Let us assume T$ = for some 
$ G X sym . Since $j G PW^ m we know that e~ 10 " lA $j(A) is bounded and analytic in the 
lower half plane, including along the rays arg A = and arg A = ir. Indeed, this follows 
from the the basic estimate 

|e-^ A $ 4 (A)| = 

for A in the lower half plane. Using the definition of T, we must have $ = — K$. The 
ith component of this equation reads 
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Using the symmetry relation $(A) = $*(A) we deduce 



Multiplying this by e 1<TiX and using that Zi = — Oie 10Li we find 



e- i<JiA $i(A)| = — ^ f 1 e ieiai ^-^- Te ~ iaj)x ^j(T)dT 



(9) 




By convexity we have 

< arg (e~ iai (mj + re iay - %)) < 7T - \oti - a^i 



(10) 



so the right hand side ([9]) is bounded for — 7r < arg A < |«j — aj_i| — n, or equivalently 

7r — \cti — < arg A < 7r. 

But since the left hand side of flU]) is also bounded along the ray arg A = 0, we can 
use Phragmen-Lindelof to deduce that e -1CTiA $j(A) is bounded in the upper half plane. 
Since we have already concluded that e -1 °" iA $j(A) is bounded in the lower half plane, it 
must be equal to a constant. This constant must be zero, however, since e 1(7iX is not 
in PW^ m . So there can be at most one element in DN($). The fact that DN($) is 
non-empty follows from Theorem 1 in [3] and its extensions. □ 

It is possible to prove that DN($) is non-empty through a direct argument by 
considering the equation T$ = A for A e Y and $ G X sym . In Lemma [6] below we 
show that T G £(X sym , Y) has closed range, so by Banach's closed range theorem it is 
sufficient to prove that A = iT* G N(J*) L for * G X sym , where T* G C(Y,X sym ) is the 
adjoint of T. For economy of presentation we leave out the straightforward argument. 

In light of the result in Lemma HI we have a well-defined map 



This is the spectral Dirichlet-Neumann map. We have the following important theorem. 

Theorem 2. The spectral Dirichlet-Neumann map DN defines a continuous linear map 
from X sym to itself. 

We will need the following simple generalisation of Lemma [2j 

Lemma 5. Any norm-bounded sequence in X sym contains a subsequence that converges 
pointwise and locally uniformly to an element ofX sym which obeys the same norm bound. 

Proof. Given a norm bounded sequence {$ m } m >i in X sym , one first applies the the result 
of Lemma[2]to the first component of the sequence to get a subsequence {$m fc }fc>i whose 
first component has the desired property. With this subsequence, one then chooses a 
sub-subsequence for which the second component has the desired properties. Continuing 
inductively gives the required result. □ 



DN : X sym > X sym : $ ^ DN($). 
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Lemma 6. The map T G £(X sym , Y) is bounded below, i.e. ||T$||y > ||$| 



x ■ 



Proof. We know that T G £(X sym , Y) and that T is continuous with respect to the 
topology of pointwise convergence. By a simple generalisation of Theorem [1] we find 
that T has closed range. Paired with the fact that T is injective on X sym , it follow^ 
that ||T$|| y > \\®\\x for $ G X sym . □ 

Proof of Theorem Linearity is straightforward so we focus on continuity. Let us 
proceed by contradiction. Suppose that the map is not continuous, i.e. it is unbounded. 
Then there exists a sequence {$ m }m>i m X sym such that ||$ m ||x = 1 for each m but 
||DN($ m )|| x — > oo. By definition we have 

T[DN($ m )] =iT[$ m ] 

for each m > 1. Introduce the new sequences 

DN($ m ) _ $ m 



l|DN($ m )|| x ' ||DN($ m )|| x 

so H^mllx = 1 f° r each m and ||A m ||x — > 0. By LemmaHwe know that T : X sym — > Y 
is continuous so we deduce ||T\I/ TO ||y — > 0. Lemma [6] gives 

1 = ||\T/ II v < HTvR \\v — )• 

which provides us with our contradiction. So there is some constant such that 
||DN($)|U< ||$|| x forall$GX sym . □ 

By multiplying the global relation by the imaginary unit i, essentially reversing the 
roles of the real and imaginary parts of X, we obtain the following. 

Corollary 1. The map DN : X sym — > X sym is a homeomorphism. 



6. Towards a New Numerical Approach 

Here we give a brief outline of how the previous theoretical results can be used to 
produce a new approach to the numerical study of the solutions to boundary value 
problems associated with Laplace's equation on the interior of a convex polygon. In 
particular, we demonstrate that this approach provides a whole family of potential 
numerical schemes that can be used to solved the global relation to any required degree 
of accuracy. A more detailed numerical study will be pursued elsewhere. 

The natural starting point is to rephrase the global relation in terms of a weak 
variational problem. Given G X sym , is it enough to find some $ G X sym such that 

T($ -i$*)(A) = 0, XeD 

\ Since R(T) C Y is closed it is also a Banach space, so T : X synl — > R(T) is a bijection between 
Banach spaces. The relevant estimate now follows from Banach's bounded inverse theorem. 
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where D is a subset of C which contains an accumulation point. This follows from a 
simple analytic continuation argument. More generally, we could study the equations 

T($-i$*)i(A) = 0, AG A, z = l,...,n (11) 

where each C C, 1 < i < n, contains an accumulation point. Motivated by these 
observations we seek to minimize the functional 



n „ 

/ |T($-in(A)| 2 ds(A) 



where 7$, 1 < i < n, are curves in C and ds(A) is the natural Lebesgue measure of arc- 
length on these curves. We impose that the curves ji are locally finite and semi-infinite, 
and each eventually coincides with the negative real axis. If we define the bilinear and 
linear forms on X sym by 



n „ 

a($, M>) = Re ]T / (T$MA)(T^MA)ds(A), 

i=l T« 
n „ 

£(*) = -lm J2 / (T$ t ),(A)(TvI/) l (A)d S (A), 

■/'Yi 



i=l 

a standard calculus of variations argument leads us to the following weak form of (II ip. 
Lemma 7. $ G X sym is a minimizer for I[<&] i/ and on/y i/ 

a($, V$el sjm . (12) 

Note that the linearity of £ and bilinearity of a follow from the fact that X sym is a real 
vector space. We need the following results to apply the standard machinery. 

Lemma 8. The bilinear form a : X sym x X sym is bounded and coercive 

(i) io($,*)i < ii^iuii^iu, (h) o($,$) > wm 



andieX: ym , i.e. < 



x ■ 



Proof. To show that a is bounded we first apply Cauchy-Schwarz 

|a(*,*)l<£ [ IO~$)i(A)||(T*)i(A)|ds(A) 

»=i 

™ / r \ 1/2 / /■ \ 1/2 

^ £ (jf l(T$) l (A)| 2 d S (A)J |(TvI/) l (A)| 2 d,(A)J 
It is now enough to show that for $ G X sym 

/ |(T$),(A)| 2 d S (A)< \\n 2 x , l<i<n. 
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Choose some R > sufficiently large so that outside the ball B R = {AgC:|A|<_R} 
all of the contours {7i}" =1 coincide with the negative real axis. We deal with the 
contributions from A G 7« D Br and A G (—00, —R) separately. We have 

2 



l(T*)i(A)p 



^i(A) iai{mi ~ m ^j(e- iA » X) 



< |$i(A)| 2 + \ eie ""'^-^^(e-^A) . 

The supremum of these terms on % D -B_r can be estimated using the standard Paley- 
Wiener inequality. Using the fact that the length of each 7$ contained in this region is 
finite, we arrive at 



(T$) 4 (A)| 2 d S (A) < 7 \\n* x . 



The contribution from (—00, —R) is easily estimated in terms of \\$\\x by using the fact 
T G C(X sym ,Y), so the claim in (i) is proven. For coercivity we note that if \\<&\\x = 1 
then it must be the case that 



n „ 

J2 / |(T$MA)| 2 d S (A)> 7 l. 

t=l T> 



(13) 



Indeed, if this were not true then there would one could construct a sequence {$ m } m >i 
with ||$ m ||x = 1 such that 



n . 1 
£ / l(T$ m ),(A)| 2 d S (A)<-. 



=1 •'T. 

Passing to a subsequence if necessary, we have $ m — >■ $' locally uniformly for some 
$' G X sym for which a($', $') = 0, i.e. (T$')j(A) = for A G 7;. By analytic continuation 
it follows that (T$')(A) = for A G C, hence $' = by the injectivity of T on X sym . We 
deduce that the sequence {$ m } m >i converges to zero locally uniformly. In particular, 
for the fixed R > used earlier and for any given e > we can take m sufficiently large 
so that 



E 

i=l 



jinB R 



-R 



(T$ m ) l (A)| 2 ds(A) 



< e. 



So for m sufficiently large we have the estimate 

a($ m ,$ m ) > ||T$ m |£-e. 

By Lemma El there is some c > such that ||T$||y > \/2c||$||_x for all $ G X sym . 
Setting e = c and choosing m sufficiently large we find 

a($ m ,$ m ) > a 
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This contradicts our assumption that a($ m , $ TO ) — > 0, so the estimate in (fT3l) must hold. 
Coercivity (ii) follows directly from f fT3|) 

n „ 

= / \(m/\m\xMX)\ 2 ds(X) > 7 ||$|&. 

i=l J H 

That £ defines a bounded linear map on X sym follows from arguments similar to those 
used to prove estimate (i). □ 

An application of the Lax-Milgram lemma gives. 

Theorem 3. There is a unique solution in X sym to the weak problem ( TJJ| ). 

Remark 5. The numerical implementation of the Fokas method has not, to date, used 
a weak approach. The standard approaches, e.g. [HI [151 [26], have approximated 
the unknown boundary boundary data {<y9 n }" =1 using suitable basis functions {$ m }m=i 
so that 

N 

<A n ( r ) ~ C *«^™( r ) l<i<n 

m=l 

for some N ^> 1. This is used to approximate the unknown spectral function $ n (A). 
Using this approximation one evaluates the global relation ([5]) at a sequence of points in 
the complex plane to get linear problem for the unknown coefficients {c.j m }. However, 
these results have been formal in nature - to this authors knowledge no proofs of 
convergence or stability have been given. 

To make rigorous these pointwise approaches it seems the semi-norm estimate 

||T$||d,oo ||$||d,oo, $eX sym , (14) 

for each open D C C is most relevant, where ||$||d,oo = maXjSup D |$j(A)|. To see this 
estimate first note the following 

||T$||d,oo >d II^IId.oo, $ e X sym and ||$|| x = 1. 

Indeed, if it were not true then one could take a sequence {$ m } m >i with ||$ m ||_D i0 o = 
II || a - = 1 an d ||T$ m ||£i )00 — > 0. Passing to a subsequence if necessary, we find 
$ m — > $ locally uniformly with T$ = on D and on all of C by analytic continuation. 
Again using the injectivity of T on X sym we deduce that $ m — >■ locally uniformly, 
contradicting our assumption that ||$ m ||D,oo — 1- So for any $ G X sym we have 

||T$||d,oc = ||$||x||T($/||$|| 

x)||d,oo ||$|| x ||($/||$|| x )|| Aoo = ||$|| A oo, 

which is the estimate in (tHJ). We note that the open set D could be replaced with any 
set containing an accumulation point. 

Using the result of Theorem [3] we can now apply standard Galerkin techniques 
to solve a sequence of finite dimensional problems whose solution approximates the 
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true solution to (1121) . the error in which is controlled by Cea's lemma [7j. A practical 
implementation of this can be achieved as follows. We write 



n 



9(\) = ^e j * j {\), 

3=1 

where {e-j}™ =1 are the usual basis vectors on R n . Since $j G PW aj for 1 < j < n, we 
may approximate each by projecting onto the finite dimensional subspace consisting of 
the sample frequencies < N. So we write 

$ ;(a) « £ $ M( A )' 

\J\<N 

where e 3 j(\) is the Jth basis function for PW aj , given explicitly by 

sin(o"jA — Ti J) 



4(A) 



CTjA — IX J 



The finite dimensional problems that approximate f lT2|) are then 



71 



/ / $^a(ej <g> ej,ej ® e}) = £(ej ® e}), 1 < i < n, \I\ < N. 

3=1 \J\<N 

This constitutes a (2N + l)n by (2iV + l)n linear system for the (2N + l)n complex 
unknowns however we are yet to take into account that $j G PWsym- We must 
enforce that $j(A) = $*(A), so that is restricted to the real part of PW aj . A 
straightforward computation reveals that this is equivalent to 

cf>- J = 

3 3 • 

By writing = Xj 7 + iY- 1 for real numbers Xj 7 and Yj 1 and noting that the previous 
symmetry condition implies = 0, we are left with a (2N + l)n by (2 AT + l)n linear 
system for the (2N + l)n real unknowns 

X°, X/, Yf, l<j<n, 1 < J < N. 

The coefficients <3>^ are then built up from the facts 

Yf J = -Yf and Ay' A'/. 



7. Conclusion 



We have shown that the global relation for Laplace's equation in a convex polygon gives 
rise to a well-defined spectral Dirichlet-Neumann map. This map is democratic, in the 
sense that the spectral boundary data $* and $ n are treated at the same level (upto 
multiplication by i). The spectral Dirichlet-Neumann map describes a homeomorphism 
on AT sym . 

We treated the case in which the Dirichlet boundary data belonged to H 1 (dQ), 
but lower regularity can be assumed. If one deals with the more general Paley- Wiener 
spaces, consisting of entire functions of exponential type whose restriction to the real 
axis belongs to L P (R) for some p > 1, then more general results can be obtained. The 
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relevant Paley-Wiener inequality in this case becomes |/(^o)| ^ eCT ' 20 ' ll/llp- It is most 
likely that these results will carry through into the limiting case where the Dirichlet 
boundary data belongs to L 2 (dQ). One would expect that the spectral Dirichlet- 
Neumann map can be extended to 

DN : FH^i-o-^G^ x • • • x JFH s+1 ^[-a ni a n ] 

^JFH s - l l 2 [-a 1 ,a 1 ]x---xFH s - l l 2 [-a n ,a n ], se [-±, §] 

so as to match the classical results for the Steklov-Poincare operator on Lipschitz 
domains [19]. Indeed, this is hinted at in the proof to Lemma H]- uniqueness followed 
from the fact that e lA<Ti = J r [5_ (7 J(A) does not belong to PW Ul = TL 2 \—<j{,ah i.e. 
S- ai L 2 = H°. However, <5_ CTi 6 H~ s for s > 1/2. These issues are discussed in [I]. 

The weak formulations given in £0 give an infinite family of problems that can be 
attacked with standard numerical procedures. Each of these problems depends on a 
choice of contours {7«}™ =1 . The constants involved in the boundedness and coercivity of 
the relevant bilinear form will depend on this choice of contours, so it is natural to ask 
which choice of contours is best. This and other aspects of the numerical implementation 
of our results are a work in progress. 

We also provided a means, via (JHJ), to make the existing numerical implementations 
of the Fokas method mathematically rigorous. 

The extension of these results to other constant coefficient elliptic boundary value 
problems is possible with suitable adjustments, but this will be pursued elsewhere. We 
expect that similar results will hold for the Helmholtz and modified Helmholtz equations, 
with the proofs following in a similar fashion. We also expect similar results to hold for 
the more important exterior problems. For an indication of the necessary modifications, 
we refer the reader to [I] where some of these modifications are presented. 

Perhaps the most important thing to note is the possible extension of the methods 
produced here to higher dimensions. It was shown in [3] that the global relation 
characterises the generalised Dirichlet-Neumann map for linear elliptic PDEs in convex 
domains in any number of dimensions. However, doing meaningful analysis with the 
global relation in higher dimensions has proved difficult, with little progress made 
over the last fifteen years. The main arguments presented here can be carried over 
to the higher dimensional problems. In particular, they can be used for boundary value 
problems in three dimensions. That this is possible is closely related to the fact that 
the theorems of Montel, Paley and Wiener (-Schwartz) extend to the complex analysis 
of several variables. 
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